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Abstract. Charged-particle production was studied in proton-proton collisions collected at the LHC with 
the ALICE detector at centre-of-mass energies 0.9 TeV and 2.36 TeV in the pseudorapidity range \rj\ < 1.4. 
In the central region (\rj\ < 0.5), at 0.9 TeV, we measure charged-particle pseudorapidity density diV c h/d?7 = 
3.02±0.01(stat)t° A ol(syst.) for inelastic interactions, and dN ch /dr] = 3.58±0.01(stat)to. 1 il(syst.) for non- 

single-diffractive interactions. At 2.36 TeV, we find dN c h/dri = 3.77 ± 0.01(stat.)^Q '^(syst.) for inelastic, 

and diVch/dr? = 4.43 ± 0.01(stat)l } ' 1 2(s?/si.) for non-single-diffractive collisions. The relative increase in 

charged-particle multiplicity from the lower to higher energy is 24.7%±0.5%(stat.)±ll%{syst.) for inelastic 

and 23.7% ± 0.5%(stat.)^\%(syst.) for non-single-diffractive interactions. This increase is consistent with 
that reported by the CMS collaboration for non-single-diffractive events and larger than that found by a 
number of commonly used models. The multiplicity distribution was measured in different pseudorapidity 
intervals and studied in terms of KNO variables at both energies. The results are compared to proton- 
antiproton data and to model predictions. 
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1 Introduction 

Whenever entering a new energy regime with hadron col- 
liders, it is important to measure the global character- 
istics of the collisions. These interactions, dominated by 
soft (i.e. small-momentum-transfer) processes, are useful 
to study QCD in the non-perturbative regime, and to 
constrain phenomenological models and event generators. 
Such studies are also important for the understanding of 
backgrounds for measurements of hard and rare interac- 
tions. 

ALICE [I] has measured the pseudorapidity density of 
charged particles produced in proton-proton collisions at 
a centre-of-mass energy y/s = 900 GeV [2] with low statis- 
tics from the first collisions at the CERN Large Hadron 
Collider (LHC) [3]. Results were given for two normaliza- 
tions: 



— inelastic (INEL); this corresponds to the sum of all in- 
elastic interactions (non-diffractive ND, single-diffrac- 
tive SD, and double-diffractive DD) with the trigger 
biases corrected for each event class individually ac- 
cording to their respective estimated abundances and 
trigger efficiencies; 

— non-singlc-diffractivc (NSD); here the corrections are 
applied to non-diffractive and double-diffractive pro- 
cesses only, while removing, on average, the single- 
diffractive contribution. 

The corrections to INEL and NSD samples are based on 
previous experimental data and simulations with Monte 
Carlo event generators. Charged-particle pseudorapidity 
density in pp collisions at LHC was also published by the 
CMS collaboration for NSD interactions [3], and by the 
ATLAS collaboration for a different event selection [5] , not 
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directly comparable with our measurements and those of 
CMS. 

Wc have used the first high energy proton-proton colli- 
sions at the LHC at a centre-of-mass energy sfs = 
2.36 TeV, as well as a larger statistics data sample at 
y/s = 0.9 TeV, to determine the pseudorapidity density of 
charged-primary particle^], dN^/dr], in the central pseu- 
dorapidity region (|^| < 1.4). According to commonly 
used models [SHU], an increase in dN c h/dr] of 17-22% for 
INEL events and of 14-19 % for NSD events is expected 
in 2.36 TeV collisions relative to 0.9 TeV collisions. 

We also studied the distribution of the multiplicity 
of charged particles in the central pseudorapidity region 
(|r;| < 1.3). The multiplicity distribution of charged par- 
ticles (the probability P(N C \ 1 ) that a collision has multi- 
plicity N c h) can be described by KNO scaling [T3] over a 
wide energy range. KNO scaling means that the distribu- 
tion (N c h)P(z), where z — N^/ (N^), is independent of 
energy. In full phase space, scaling holds up to the top ISR 
energy (pp at y/s = 62.2 GeV) [14]. Deviations from scal- 
ing are observed at higher energies, starting at 200 GeV 
with pp collisions at the SppS collider [15] . However, in 
limited central 77-intervals scaling has been found to hold 
up to 900 GeV. The UA5 collaboration [TB] observed scal- 
ing for non-single-diffractive events in restricted central 
^-intervals and its progressive violation with increasing rj- 
ranges. The UA1 collaboration [T7] also observed scaling 
in a larger interval \r]\ < 2.5. In inelastic events, devia- 
tion from KNO scaling was observed in full phase space 
already at ISR energies [T3] . Such deviations are generally 
attributed to semi-hard gluon radiation (minijets) and to 
multi-parton scattering. 

The Negative-Binomial Distribution (NBD) [H] de- 
scribes multiplicity distributions in full phase space up 
to 540 GeV; however, this description is not successful 
at 900 GeV [T5]. NBD describes the distributions up to 
1.8 TeV in limited ry-intervals (\t]\ < 0.5) [5D]- For larger 
//-intervals and in full phase space, only the sum of two 
NBDs provides a reasonable fit [2"T1|2"2"] . 

Comparing these multiplicity measurements with the 
predictions of Monte Carlo generators used by the LHC 
experiments will allow a better tuning of these models to 
accurately simulate minimum-bias and underlying-event 
effects. A recent review of multiplicity measurements at 
high energies can be found in [2"3"] . 

This article is organized as follows: a description of 
the ALICE detector subsystems used in this analysis is 
presented in Section 2; Section 3 is dedicated to the def- 
inition of the event samples; Section 4 to data analysis; 
in Section 5 systematic uncertainties are discussed; the 
results are given in Section 6 and Section 7 contains the 
conclusions. 



2 The ALICE experiment and data collection 

The ALICE experiment consists of a large number of de- 
tector subsystems which are described in detail in [I] . This 
analysis is based mainly on data from the Silicon Pixel 
Detector (SPD), since it has the largest pseudorapidity 
coverage in the central region and is located closest to the 
interaction region, implying a very low momentum cut-off 
and a small contamination from secondary particles. 

The SPD detector surrounds the central beryllium 
beam pipe (3 cm radius, 0.23 % of a radiation length) with 
two cylindrical layers (at radii of 3.9 and 7.6 cm, 2.3% of 
a radiation length) and covers the pseudorapidity ranges 
\rj\ < 2 and < 1.4 for the inner and outer layers, respec- 
tively. The number of inactive (dead or noisy) individual 
pixels is small, about 1.5%, but in addition some 17% of 
the total area is currently not active, mostly because of 
insufficient cooling flow in some of the detector modules. 
The number of noise hits in the active pixels of the SPD 
was measured with a random trigger to be of the order 
of 10 -4 per event. The SPD was aligned using cosmic-ray 
tracks [23] collected prior to the collider run and tracks 
from collisions recorded at ^/s — 0.9 TeV. 

Information from two scintillator hodoscopes, called 
VZERO counters, was used for event selection and back- 
ground rejection. These counters are placed on either side 
of the interaction region at z = 3.3 m and z — —0.9 m. 
They cover the regions 2.8 < r\ < 5.1 and —3.7 < 77 < —1.7 
and record both amplitude and time of signals produced 
by charged particles. 

The central detector subsystems are placed inside a 
large solenoidal magnet which provides a field of 0.5 T. For 
the 2.36 TeV data taking the VZERO detectors were not 
turned on. Therefore, the trigger conditions, the analysis 
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in the collision and all decay products, except products from 
weak decays of strange particles. 



Fig. 1. Distributions of reconstructed event vertices along the 
beam direction (z) obtained from hit correlations in the two 
pixel layers of the ALICE inner tracking system for the event 
samples used in the analysis (see text): ^J~s — 0.9 TeV (full 
symbols) and yfs = 2.36 TeV (open symbols). The lines are 
from Monte Carlo simulations. Vertical dashed lines delimit the 
region \z\ < 10 cm, where the events for the present analysis 
were selected. 
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and the systematic errors differ slightly between the two 
data sets (see below). 

Because of the low interaction rate it was possible to 
use a rather loose trigger for collecting data. At 0.9 TeV, 
the minimum-bias trigger required a hit in either one of 
the VZERO counters or in the SPD detector; i.e. essen- 
tially at least one charged particle anywhere in the 8 units 
of pseudorapidity covered by these trigger detectors. At 
the higher energy, the trigger required at least one hit in 
the SPD detector (|?7| < 2). The events were collected in 
coincidence with the signals from two beam pick-up coun- 
ters (BPTX), one on each side of the interaction region, 
indicating the presence of passing bunches. 

The bunch intensity was typically 5 x 10 9 protons, giv- 
ing a luminosity of the order of 10 26 cm _2 s^ 1 . This value 
corresponds to a rate of a few Hz for inelastic proton- 
proton collisions and a negligible pile-up probability for 
events in the same bunch crossing. 

In the case of the 0.9 TeV data, events in coincidence 
with only one passing bunch, as well as when no bunch was 
passing through the detector, were also registered. These 
control triggers were used to measure the beam-induced 
and accidental backgrounds. 

The observed longitudinal sizes of the interaction re- 
gions can be inferred from Fig. [TJ which shows the dis- 
tribution of the interaction vertices along the beam axis 
reconstructed from hit correlation in the two SPD lay- 
ers. The vertex distribution has an estimated r.m.s. of 
4.1cm and 2.7cm for the 0.9 TeV and 2.36 TeV samples, 
respectively. These vertex distributions are for all trig- 
gered events with a reconstructed primary vertex after 
background rejection. They are compared to Monte Carlo 
simulations of proton-proton collisions using a Gaussian 
beam profile with a standard deviation as measured from 
data. The experimentally observed tails (mainly from 
events with one or two reconstructed tracks) are well de- 
scribed in the simulation, confirming that beam induced 
background is very small in the selected sample. The ver- 
tex-reconstruction efficiency is practically independent of 
the vertex z-position for \z\ < 10 cm. 



3 Event selection and corrections to IN EL 
and NSD event classes 

Slightly different event selections were applied after data 
reconstruction for the analysis of the two collision energies 
because of the different detector configurations. 

For both data samples, an offline selection is applied to 
reject beam induced background. At 0.9 TeV, the VZERO 
counters were used to remove beam-gas or beam-halo 
events by requiring their timing signals, if present, to be 
compatible with particles produced in collision events (see 
[2] for more details). At both energies, this background 
was also rejected by exploiting the correlation between the 
number of clusters of pixel hits and the number of track- 
lets (short track segments in the SPD, compatible with 
the event vertex, as described below). From the analysis 
of our control triggers, we found that background events 
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Fig. 2. Charged-particle multiplicity dependence of the se- 
lection efficiency, for INEL events (triangles) with MBor se- 
lection, and for NSD (circles) and SD (squares) events with 
MBand selection obtained with a simulation of the 0.9 TeV 
data. 



typically have a large number of pixel hits compared with 
the number of tracklets pointing to the reconstructed ver- 
tex. 

At 0.9 TeV, for the INEL analysis, we used the trig- 
gered event sample requiring a logical OR between the sig- 
nals from the SPD and VZERO detectors (MB r). How- 
ever, for the NSD analysis we selected a subset of the total 
sample by requiring a coincidence between the two sides of 
the VZERO detectors (MBand)- This requires the detec- 
tion of at least one charged particle in both the forward 
and backward hemispheres, which are separated by 4.5 
units of pseudorapidity. In this subset, single-diffraction 
events are suppressed, therefore, model dependent correc- 
tions and associated systematic errors are reduced (see 
below). The selection efficiencies, MBor f° r INEL events 
and MBand for NSD events, are multiplicity dependent 
as illustrated in Fig. [2] As expected, the MBand selec- 
tion has a low efficiency for SD events, in particular at 
low multiplicities, where they contribute most. After these 
selections, the remaining background at 0.9 TeV was es- 
timated, and corrected for, with the help of the control 
triggers. The background events (99 % of which have no 
reconstructed tracklets) correspond to about 2 % of the 
events in the INEL sample and to less than 0.01 % in the 
NSD sample. 

The 2.36 TeV data sample was triggered by at least one 
hit in the SPD (MBspd) and this selection was used for 
both INEL and NSD analyses. After rejecting the back- 
ground using the correlation between the number of pixel 
hits and the number of tracklets, the remaining back- 
ground (93 % of which has no reconstructed tracklets) was 
estimated to be 0.7%. We have assumed that the corre- 
lation between the number of clusters of pixels and the 
number of tracklets is similar at both energies because 
accelerator and detector conditions did not change signif- 
icantly between the two data collection periods. 
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Table 1. (a) Relative fractions of SD and DD events, as obtained from previous measurements at 0.9 TeV [27] and 1.8 
TeV [301133] . The measured DD fractions are scaled according to the prescription in [33] . Corresponding fractions calculated 
using PYTHIA and PHOJET are given for events within the diffractive-mass range covered experimentally (see text), and also 
without the restriction on diffractive-mass (parentheses), (b) Selection efficiencies for different classes of events: at 0.9 TeV, 
where the MBor selection was used for INEL sample and MB and for NSD sample; at 2.36 TeV, where the selection using the 
SPD only was used for both INEL and NSD samples. 
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SD 
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0.159 ±0.024 


0.167 (0.209) 


0.126 (0.161) 


DD 


0.095 ± 0.060 
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0.066 
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(b) Selection 


efficiencies 
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2.36 TeV 




PYTHIA 


PHOJET 


PYTHIA 


PHOJET 




MBor 


MB AND 


MBor 


MB AND 


MBspd 


MBspd 


SD 


0.77 


0.29 


0.86 


0.34 


0.55 


0.62 


DD 


0.92 


0.49 


0.98 


0.77 


0.63 


0.79 


ND 


1.00 


0.98 


1.00 


0.96 


0.99 


0.99 


INEL 


0.95 




0.97 




0.86 


0.90 


NSD 




0.92 




0.94 


0.94 


0.97 



In both data samples, the cosmic-ray contamination, 
estimated from the control triggers and from absolute 
rates, is negligible. Additional crosschecks of background 
levels were made by visual scanning of a few hundred se- 
lected events. 

The number of collision events used in this analysis 
corresponds to about 150 000 and 40 000 interactions for 
the 0.9 and 2.36 TeV data, respectively. 

The efficiencies of our selections and their sensitivi- 
ties to variations in the relative fractions of event classes 
were studied using two different Monte Carlo generators, 
PYTHIA 6.4.14 and 6.4.21 [7J|B] tune D6T [9] and PHO- 
JET 1.12 [12] used with PYTHIA 6.2.14, with the detec- 
tor simulation and reconstruction software framework Ali- 
Root [35], which includes a detailed model of the ALICE 
apparatus. Particle transport in the detector was simu- 
lated using the GEANT-3 software package [2"5] . 

To normalize our results to INEL and NSD event class- 
es, we used measured cross sections in pp collisions instead 
of those provided by the event generators. The single- 
diffraction production cross section at 0.9 TeV was mea- 
sured by UA5 [37] in the kinematic range M 2 < 0.05s, 
where M is the mass of the diffractive system. The pub- 
lished value is a result of an extrapolation to low diffractive 
masses M < 2.5 GeV and therefore is model-dependent. 
There exist indications from measurements at lower en- 
ergy, y/s = 546 GeV [25], and from phenomenological 
models [2j|] that this cross section value may be under- 
estimated by up to 30 %. We decided to use the published 
value, but we checked that our results stay within their 
systematic uncertainties, if instead we use a 30% higher 
single-diffractive cross section. 

At 1.8 TeV, the single-diffraction cross section was 
measured by E710 [30] and CDF [31] in the kinematic 
ranges 2 GeV 2 < M 2 < 0.05s and 1.4 GeV 2 < M 2 < 
0.15s, respectively. As the CDF result includes significant 



model dependent acceptance corrections at low masses, we 
used the E710 measurement. The inelastic cross section at 
1.8 TeV, needed for normalization, was taken from |32) . 
At both energies, the fraction of SD events in the Monte 
Carlo generators was normalized to the data in the mass 
regions covered by the corresponding experiments. 

To simulate the SD events we used the two Monte 
Carlo event generators and rely on them to tag these 
events. In order to classify an event as diffractive, the 
diffractive mass should satisfy a coherence condition (M 2 
< Ks, where K = 0.05-0.15), which effectively provides a 
high-mass cut-off. PHOJET uses this coherence condition 
with K — 0.15 and we checked that, if we would further 
restrict diffractive masses using K = 0.05, our multiplic- 
ity results will practically not change. PYTHIA generates, 
with a low probability, SD events with very high diffractive 
masses. Therefore, we checked the stability of our results 
by imposing the tight coherence condition (K — 0.05) to 
generated SD events. This decreased the average multi- 
plicities for the NSD samples by less then 2 %, well within 
our systematic uncertainties. 

Measurements of double-diffraction cross sections are 
available from UA5 [27] at 0.9 TeV and CDF [33] at 
1.8 TeV. Experimentally, DD events are defined by re- 
quiring a minimum pseudorapidity gap (of about 3 units), 
where no charged particles are detected. When implement- 
ing these experimental cuts in the event generators, the 
results were widely fluctuating and inconsistent with the 
measurements, possibly because the occurrence of large 
rapidity gaps is very sensitive to the model assumptions 
and process parameterizations. Therefore, for classifica- 
tion of DD events we used the process type information 
provided by the generators but we adjusted the fractions 
to the measured values. The values used take into account 
an increase of the DD fractions due to the pseudorapidity- 
gap definition as described in [33] . Note that the correction 
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arising from unmeasured DD events is small, both because 
the cross section for DD is small and because the event 
selection efficiency is large in our samples. 

The relative fractions of SD and DD events, as mea- 
sured in 27,30,33 arc summarized in Table [TJ along with 
our calculated trigger and selection efficiencies. The rela- 
tive fractions for SD and DD vary very slowly with energy, 
therefore, we used the measurements available at 1.8 TeV 
for the 2.36 TeV sample. 



4 Analysis method 

The analysis method is based on using hits in the two 
SPD layers to form short track segments, or tracklets. This 
method is similar to that used by the PHOBOS experi- 
ment with the first heavy-ion data from RHIC [33]. We 
start with the reconstruction of the position of the inter- 
action vertex by correlating hits in the two silicon-pixel 
layers. The vertex resolution achieved with this simple 
method depends on the track multiplicity, and is typically 
0.1-0.3 mm in the longitudinal (z) and 0.2-0.5 mm in the 
transverse direction. For events with only one SPD track- 
let, the z- vertex position is determined by the point of 
closest approach to the mean beam axis. A vertex was re- 
constructed for 83 % of events in the MBor selection and 
for 93 % of events in the MB and selection. At the higher 
energy, in the MBspd selection 93 % of events have a ver- 
tex reconstructed. Events with vertices within \z\ < 10 cm 
are used in this analysis. 

Using the reconstructed vertex as the origin, we calcu- 
late the differences in azimuthal (Aip, bending plane) and 
polar (A9, non-bending direction) angles of two hits, one 
in the inner and one in the outer SPD layer. Hit combina- 
tions, called tracklets, are selected by a cut on the sum of 
the squares of Aip and A9, with a cut-off at 80mrad and 
25 mrad, respectively. The cut imposed on the difference 
in azimuthal angles would reject charged particles with 
a transverse momentum below 30MeV/c; however, the 
effective transverse-momentum cut-off is determined by 
particle absorption in the material and is approximately 
50 MeV/c. If more than one hit in a layer matches a hit in 
the other layer, only the hit combination with the smallest 
angular difference is used. 

For the pseudorapidity-density measurement, all events 
with vertex in the range \z\ < 10 cm are used. For multi- 
plicity-distribution measurements, the whole 77-interval 
considered has to be covered by the acceptance of the 
SPD, for every event. Therefore, only events from a limited 
z-range of collision vertices are used for the two largest rj- 
intervals, which reduces the available event statistics. At 
0.9 TeV these reductions are 15% for \rj\ < 1.0 and 60% 
for \r)\ < 1.3, and at 2.36 TeV 4% for \r)\ < 1.0 and 46% 
for |ry| < 1.3. 

The number of primary charged particles is estimated 
by counting the number of tracklets. This number was 
corrected for: 

— geometrical acceptance, detector and reconstruction 
efficiencies; 



— contamination by weak-decay products of long-lived 
particles (K°, A, etc.), gamma conversions and sec- 
ondary interactions; 

— undetected particles below the 50MeV/c transverse- 
momentum cut-off; 

— combinatorial background caused by an accidental as- 
sociation of hits in the two SPD layers, estimated from 
data by counting pairs of hits with a large Aip. 

The probability of an additional collision in the same 
bunch crossing (pile-up) at the estimated luminosity is 
below 10 -3 . The effect on both multiplicity density and 
multiplicity distribution measurements due to such events 
has been found to be negligible. Particular attention was 
paid to events having zero or one charged tracklets in the 
SPD acceptance. For the 0.9 TeV sample, the number of 
zero-track events for \z\ < 10 cm was estimated using trig- 
gered events without a reconstructed vertex. At 2.36 TeV, 
due to the different trigger (see Section 2), we have to 
use Monte Carlo simulations to estimate this number and 
therefore the results are more model-dependent than those 
at 0.9 TeV. As a consequence, the size of systematic un- 
certainties on average multiplicity is bigger at 2.36 TeV 
than that at 0.9 TeV, as described in Section 5. 

The total number of collisions used for the normaliza- 
tion was calculated from the number of events with recon- 
structed vertex selected for the analysis and the number 
of triggered events without vertex. The latter number was 
corrected for beam-induced and accidental background as 
measured with the control triggers (see Section 2). A small 
correction, determined from simulations, is applied to the 
number of tracks due to events with no reconstructed ver- 
tex. In order to get the normalization for INEL and NSD 
events, we further corrected the number of events for the 
selection efficiency for these two event classes. For NSD 
events, we subtracted the single-diffractive contribution. 
The selection efficiencies depend on the charged-particle 




Measured multiplicity N, 



Fig. 3. Measured raw multiplicity distributions in three pseu- 
dorapidity ranges for both energy samples (0.9 TeV full lines, 
2.36 TeV dashed lines). Note that for \t)\ < 1.0 and |r?| < 1.3 
the distributions have been scaled for clarity by the factor in- 
dicated. 
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Fig. 4. Graphical representation of the detector response 
matrix: number of tracklets found in the SPD (m) vs the 
number of generated primary particles in 77 j < 1.0 (t) for 
^/s — 0.9 TeV. The distribution of the measured tracklet mul- 
tiplicity for a given generated multiplicity shown with its most 
probable value (dots), r.m.s. (shaded areas), and full spread 
(squares) . 



multiplicity, as shown in Fig. [2] for the 0.9 TeV data sam- 
ple for different event classes (INEL, NSD, and SD). At 
both energies, the efficiency is close to 100% for multi- 
plicities of one or above for the INEL class, and reaches 
90 % for multiplicities above two for the NSD class. The 
averaged combined corrections in number of events due to 
the vertex-reconstruction and the selection efficiencies for 
INEL collisions are 5 % and 24 % for 0.9 TeV and 2.36 TeV 
data, respectively. This correction is larger at the higher 
energy because of significantly smaller pseudorapidity cov- 
erage of the MBspd selection compared with the MBor 
selection and the necessity for large correction for zero- 
multiplicity events at this energy. For NSD collisions, at 
both energies, these event-number corrections are small 
(2% and 1% for 0.9 TeV and 2.36 TeV data, respec- 
tively) as a consequence of partial cancelation between 
adding non-observed ND and DD events, and subtracting 
triggered SD events. The resulting model-dependent cor- 
rection factors due to the selection efficiencies applied to 
averaged charged-particle multiplicities for the NSD sam- 
ples are 0.973 and 1.014 for 0.9 TeV and 2.36 TeV data, 
respectively. 

The multiplicity distributions, measured in three r\- 
intervals, are shown in Fig. [3] for raw data at both ener- 
gies. The method used to correct the raw measured dis- 
tributions for efficiency, acceptance, and other detector 
effects, is based on unfolding with \ 2 minimization with 
regularization [35] . The detector response was determined 
with the same Monte Carlo simulation as described above. 
Figure 0] illustrates the detector response matrix R mt for 
1 77 1 < 1, which gives the conditional probability that a 
collision with multiplicity t is measured as an event with 
multiplicity m. Therefore, each column is normalized to 
unity. This matrix characterizes the properties of the de- 
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Fig. 5. Measured raw multiplicity distribution (elements of 
vector M, histogram), superimposed on the convolution RU of 
the unfolded distribution with the response matrix (crosses), 
at \fs = 0.9 TeV for < 1.0 (upper plot). The error bars are 
omitted for visibility. Normalized residuals, i.e. the difference 
between the measured raw distribution and the corrected dis- 
tribution folded with the response matrix divided by the mea- 
surement error (lower plot). The inset shows the distribution of 
these normalized residuals fitted with a Gaussian (ex ~ 1.06). 



tector and does not depend on the specific event gener- 
ator used for its determination, apart from second-order 
effects due to, for example, differences in particle com- 
position and momentum spectra, discussed in Section 5. 
As this matrix is practically independent of energy, it is 
shown for the 0.9 TeV case only. The unfolded spectrum 
U{N c \ l ) is found by minimizing 



x 2 (t/) = E 



+ PF(U), (1) 



where R is the response matrix, M is the measured spec- 
trum, e is the estimated measurement error, and f3F(U) is 
a regularization term that suppresses high-frequency com- 
ponents in the solution. The only assumption made about 
the shape of the corrected spectrum is that it is smooth. 
The smoothness is imposed by the choice 



Ml! 

Ut 



E 



(Ut-i - Ut) 2 
Ut 



(2) 



which minimizes the fluctuations with respect to a con- 
stant constraint imposed by first derivatives. The regular- 
ization coefficient /3 is chosen such that, after minimiza- 
tion, the contribution of the first term in Eq. [1] is of the 
same order as the number of degrees of freedom (the num- 
ber of bins in the unfolding). 

The unfolded spectrum is corrected further for vertex 
reconstruction and event selection efficiencies (see Fig. [5]) 

The behaviour of the deconvolution method is illus- 
trated in Fig. [5] for the case \r)\ < 1.0 at y/s = 0.9 TeV 
showing that the normalized residuals are well-behaved 
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Table 2. Contributions to systematic uncertainties in the measurements of the charged-particle pseudorapidity density and of 
the multiplicity distribution. For pseudorapidity densities, when two values are given, they correspond to the pseudorapidities 
0.0 and 1.4, respectively. The sign of the event-generator uncertainties indicates if the result using PHOJET corrections is higher 
(positive sign) or lower (negative sign) than that using PYTHIA corrections. For multiplicity distributions the values are given 
for \rj\ < 1.0. Multiple values indicate uncertainties for respective multiplicities shown in parentheses. 



Uncertainty 


diVoh/d?? analysis 


P(N c h) analysis 






0.9 TeV 


2.36 TeV 


0.9 TeV 


2.36 TeV 


Tracklet selection cuts 


negl. 


negl. 


negl. 


negl. 


Material budget 


negl. 


negl. 


negl. 


negl. 


Misalignment 


negl. 


negl. 


negl. 


negl. 


Particle composition 


0.5-1.0 % 


0.5-1.0 % 


included in detector efficiency 


Transverse-momentum spectrum 


0.5% 


0.5% 


included in detector efficiency 


Contribution of diffraction (INEL) 


0.7% 


2.6% 


3-0 % (0-5) 


5-0 % (0-5) 


Contribution of diffraction (NSD) 


2.8% 


2.1% 


24-0 % (0-10) 


12-0% (0-10) 


Event-generator dependence (INEL) 


+1.7% 


+5.9 % 


8-0 % (0-5) 


25-0 % (0-10) 


Event-generator dependence (NSD) 


-0.5% 


+2.6 % 


3-5-1 % (0-10-40) 


32-8-2 % (0-10-40) 


Detector efficiency 


1.5% 


1.5% 


2-4-15 % (0-20-40) 


3-0-9 % (0-8-40) 


SPD triggering efficiency 


negl. 


negl. 


negl. 


negl. 


VZERO triggering efficiency (INEL) 


negl. 


n/a 


negl. 


n/a 


VZERO triggering efficiency (NSD) 


0.5% 


n/a 


1% 


n/a 


Background events 


negl. 


negl. 


negl. 


negl. 


Total (INEL) 


+2.5 0/ 
-1.8 /o 


+6.7 0/ 
-3.1 /o 


9-4-15 % (0-20-40) 


25-0-9 % (0-10-40) 


Total (NSD) 


+3.3 <w 
-3.3 /o 


+3.7 0/ 
-2.7 S° 


24-5-15 % (0-10-40) 


32-8-9 % (0-10-40) 



over the whole measured multiplicity range. The \ 2 dif- 
ference between the measured raw distribution and the 
corrected distribution folded with the response matrix is 
X 2 /ndf = 36.7/35 = 1.05. Similar behaviour is observed 
for other rj intervals and at 2.36 TeV. 

We checked the sensitivity of our results to - 

— The value of the regularization coefficient /3. 

— Changing the regularization term, defined in Eq. [21 to: 



{U't'f 



E 



(Ut-! - 2U t + U t . 



-i) 2 



Ui 



(3) 



which minimizes the fluctuations with respect to a lin- 
ear constraint imposed by second derivatives. 
Changing the unfolding procedure. An unfolding based 
on Bayes' theorem [3"6ll3"T] produces consistent results. 
It is an iterative procedure using the relations: 

Rmt • Pt 



Rfm — 



Y]fi Rmt'Pf 



u t = j2 Rt 



(4) 



with an a priori distribution P. The result U of an 
iteration is used as a new a priori P distributions for 
the following iteration. 
— Variation of convergence criteria and initial distribu- 
tion. For both unfolding procedures we checked that 
the results are insensitive to the details of the conver- 
gence criteria and a reasonable choice of initial distri- 
butions. 

The details of this analysis are described in [55] . 

5 Systematic uncertainties 

In order to estimate the systematic uncertainties, the abo- 
ve analysis was repeated: 



— varying the A<p and AO cuts used for the tracklet def- 
inition by ±20%; 

— varying the density of the material in the tracking sys- 
tem, thus changing the material budget by ±10%; 

— allowing for detector misalignment by an amount of up 
to 100 //m; 

— varying the composition of the produced particle types 
with respect to the yields suggested by the event gen- 
erators by ±30%; 

— varying the non-observed-particle yield below the trans- 
verse momentum cut-off for tracklet reconstruction by 
±30%; 

— varying the ratios of the ND, SD, and DD cross sec- 
tions according to their measured values and errors 
shown in Table [TJ thus evaluating the uncertainty in 
the normalization to INEL and NSD events; 

— varying the thresholds applied to VZERO counters, 
both in simulation and in data (for the 0.9 TeV sam- 
pie). 

The results are summarized in Table [2] using the correc- 
tions calculated with PYTHIA tune D6T. Whenever cor- 
rections obtained with PHOJET give a different value, the 
difference is used in calculating an asymmetric systematic 
uncertainty. These two models were chosen because they 
predict respectively the lowest and the highest charged- 
particle densities for INEL collisions at both energies (see 
Section 6). 

The SPD efficiencies for trigger and for pixel hits are 
determined from the data. The SPD trigger efficiency is 
determined to be 98% with negligible uncertainty based 
on analysis of the trigger information recorded in the data 
stream for events with more than one tracklet. The detec- 
tor efficiency is determined from pixel-hit distributions, 
and checked by tracklet reconstruction. The uncertainty 
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Fig. 6. Left: measured pseudorapidity dependence of AN C ^/Arj at *Js = 0.9 TeV for INEL (full symbols) and NSD (open 
symbols) collisions. The ALICE measurements (squares) are compared to UA5 pp data 40 (triangles) and to CMS pp data at 
the LHC [4j (stars). Right: measured pseudorapidity dependence of AN^/Ar] at *Js = 2.36 TeV for INEL (full symbols) and 
NSD (open symbols) collisions. The ALICE measurement (squares) for NSD collisions is compared to CMS NSD data 4 (stars) 
and to model predictions, PYTHIA tune D6T [5] (solid line) and PHOJET [TJ] (dashed line). For the ALICE data, systematic 
uncertainties are shown as shaded areas; statistical uncertainties are invisible (smaller than data marks). For CMS data error 
bars show the statistical and systematic uncertainties added in quadrature. 



on the detector acceptance and efficiency due to the lim- 
ited hit statistics and the current alignment precision of 
the detector is estimated by this method to be 1.5%. The 
uncertainty in background corrections was estimated ac- 
cording to the description in Section 3. 

The total systematic uncertainty on the pseudorapid- 
ity density measurement at 0.9 TeV is smaller than 2.5% 
for INEL collisions and is about 3.3% for NSD collisions. 
At 2.36 TeV, the corresponding uncertainties are below 
6.7% and 3.7% for INEL and NSD collisions, respectively. 
For all cases, they are dominated by uncertainties in the 
cross sections of diffractive processes and their kinematics. 

To evaluate the systematic error on the multiplicity 
distribution, a new response matrix was generated for each 
change listed above and used to unfold the measured spec- 
trum. The difference between these unfolded spectra and 
the unfolded spectrum produced with the unaltered re- 
sponse matrix determines the systematic uncertainty. 

Additional systematic uncertainties originate from the 
unfolding method itself, consisting of two contributions. 
The first one arises from statistical fluctuations due to 
the finite number of events used to produce the response 
matrix as well as the limited number of events in the 
measurement. The unfolding procedure was repeated 100 
times while randomizing the input measurement and the 
response matrix according to their respective statistical 
uncertainties. The resulting uncertainty due to the re- 
sponse matrix fluctuations is negligible. The uncertainty 
on the measured multiplicity distribution due to the event 
statistics reproduces the uncertainty obtained with the 
minimization procedure, as expected. 

A second contribution arises from the influence of the 
regularization on the distribution. The bias introduced by 
the regularization was estimated using the prescription de- 
scribed in 39] and is significantly lower than the statistical 



error inferred from the x 2 minimization, except in the low- 
multiplicity region. In that region, the bias is about 2 %, 
but the statistical uncertainty is negligible. Therefore, we 
added the estimated value of the bias to the statistical 
uncertainty in this region. The correction procedure is in- 
sensitive to the shape of the multiplicity distribution of 
the events, which produce the response matrix. 

Table [2] summarizes the systematic uncertainties for 
the multiplicity distribution measurements. Note that the 
uncertainty is a function of the multiplicity which is re- 
flected by the ranges of values. Further details about the 
analysis, corrections, and the evaluation of the systematic 
uncertainties are in |38j . 

Both the pseudorapidity density and multiplicity dis- 
tribution measurements have been cross-checked by a sec- 
ond analysis employing the Time-Projection Chamber 
(TPC) pQ. It uses tracks and vertices reconstructed in the 
TPC in the pseudorapidity region \rj\ < 0.8. The pseu- 
dorapidity density is corrected using a method similar to 
that used for the SPD analysis. The results of the two 
independent analyses are consistent. 



6 Results 

In this section, pseudorapidity density and multiplicity 
distribution results are presented for two centre-of-mass 
energies and compared to results of other experiments 
and to models. For the model comparisons we have used 
QGSM [5], three different tunes of PYTHIA, tune D6T [5], 
tune ATLAS-CSC [TDj and tune Perugia-0 [TJ, and PHO- 
JET [12]. The PYTHIA tunes have been developed by 
three independent groups extensively comparing Monte 
Carlo distributions to underlying-event and minimum-bias 
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Table 3. Charged-particle pseudorapidity densities measured by ALICE in the central pseudorapidity region (\rj\ < 0.5), for 
inelastic (INEL), non-single-diffractive (NSD), and inelastic with > (INEL>0) proton-proton collisions at centre-of- 
mass energies of 0.9 TeV and 2.36 TeV. The ratios of multiplicity densities between the two energies are also given. Data at 
yfs = 0.9 TeV are compared to CMS NSD data g] and UA5 NSD and INEL pp data gD]. Data at yfs = 2.36 TeV are compared 
to CMS NSD data. For ALICE and CMS measurements, the first error is statistical and the second one is systematic; no 
systematic uncertainty is quoted by UA5. These data are also compared to predictions for pp collisions from different models: 
QGSM _6 , PYTHIA tune D6T [5] (a), tune ATLAS-CSC 10 (6), and tune Perugia-0 [II] (c), and PHOJET [12]. 



Experiment 


ALICE pp 


CMS pp 




UA5 pp 


QGSM 




PYTHIA 




PHOJET 


Model 












a 


b 


c 










= 0.9 TeV 














INEL 


3.02 ± 0.01±o'.o5 






3.09 ±0.05 


2.98 


2.35 


3.04 


2.46 


3.21 


NSD 


3.58 ± 0.01±g i2 


3.48 ±0.02 ±0.13 


3.43 ± 0.05 


3.47 


2.85 


3.74 


3.02 


3.67 


INEL>0 


4.20 ±0.01 ± 0.03 










3.40 


4.35 


3.61 


4.06 








= 2.36 TeV 














INEL 


3.77 ± 0.0llo.?2 








3.65 


2.81 


3.64 


2.94 


3.76 


NSD 


4.43 ± O.Ollo'ia 


4.47 ±0.04 ±0.16 




4.14 


3.38 


4.44 


3.57 


4.20 


INEL>0 


5.13 ±0.03 ±0.03 










3.95 


5.05 


4.18 


4.62 


Ratios 


INEL 


1.247 ± 0.005±o.o2g 








1.22 


1.20 


1.20 


1.20 


1.17 


NSD 


1.237 ±0.005±oon 


1.28 ±0.014 ±0.026 




1.19 


1.19 


1.19 


1.18 


1.14 


INEL>0 


1.226 ±0.007 ± 0.010 










1.16 


1.16 


1.16 


1.14 



Tevatron data. Data from hadron colliders at lower en- 
ergies have been used to fix the energy scaling of the 
parameters. Tune D6T uses the old PYTHIA multiple 
scattering and Q 2 -ordered showers, whereas the two other 
tunes use the new multiple-scattering model provided by 
PYTHIA 6.4 and transverse-momentum-ordered shower- 
ing. Perugia-0 was not tuned for diffractive processes, 
which affects the validity of this tune for the lowest multi- 
plicities. For final-state-radiation and hadronization, Peru- 
gia-0 adds parameters fitted to LEP data. The charged- 
particle density in the central rapidity region is mainly 
influenced by the infrared cut-off for parton scattering at 
the reference energy (1.8 TeV) and its energy dependence. 

Figure [5] (left) shows the charged-particle density as 
a function of pseudorapidity obtained for INEL and NSD 
interactions at a centre-of-mass energy ^/s = 0.9 TeV com- 
pared to pp data from the UA5 experiment [40] . and to 
pp NSD data from the CMS experiment [4]. The result 
is consistent with our previous measurement 2 and with 
UA5 and CMS data. Figure [5] (right) shows the measure- 
ment of d-/V c h/d?7 for INEL and NSD interactions at yfs = 
2.36 TeV compared to CMS NSD data g] and to PYTHIA 
tune D6T and PHOJET calculations. Our results for NSD 
collisions are consistent with CMS measurements, system- 
atically above the PHOJET curve, and significantly higher 
than the distribution obtained with the PYTHIA tune 
D6T. Note that in the CMS pseudorapidity-density mea- 
surement the contribution from charged leptons was ex- 
cluded. This implies that the CMS value is expected to 
be approximately 1.5% lower than in our result, where 
charged leptons are counted as primary particles. 

The pseudorapidity density measurements in the cen- 
tral region (|?7| < 0.5) are summarized in Table [3] along 
with model predictions obtained with QGSM, PHOJET 
and three different PYTHIA tunes. Note that QGSM is 



not readily available as an event generator and the pre- 
dictions for some of the event classes were obtained ana- 
lytically by the authors of [6]. At both energies, PYTHIA 
tune D6T and PHOJET yield respectively the lowest and 
highest charged-particle densities for INEL collisions. 

Because part of the systematic uncertainties cancels 
in the ratio of the multiplicity densities between the two 
energies, these ratios are compared to model calculations 
as well. The main contribution to the systematic uncer- 
tainties in the measurement of charged-particle densities 




Energy \[s (GeV) 

Fig. 7. Charged-particle pseudorapidity density in the central 
rapidity region in proton-proton and proton-antiproton inter- 
actions as a function of the centre-of-mass energy. The dashed 
and solid lines (for INEL and NSD interactions, respectively) 
show a fit with a power-law dependence on energy. Note that 
data points at the same energy have been slightly shifted hor- 
izontally for visibility. 
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Fig. 8. Corrected multiplicity distributions in three pseudorapidity ranges for NSD events. The solid lines show NBD fits. Error 
bars represent statistical uncertainties and shaded area systematic ones. Left: data at yfs = 0.9 TeV. The ALICE measurement 
for \r}\ < 0.5 is compared to the UA5 data at the same energy |19| . In the inset the ratio of these two measurements is shown, 
the shaded area represents our combined statistical and systematic uncertainty, and the error bars those of UA5. Right: data 
at yfs — 2.36 TeV. Note that for \rj\ < 1.0 and \rj\ < 1.3 the distributions have been scaled for clarity by the factor indicated. 



comes from the estimate of the number of events with 
zero tracks. Therefore, in addition to the two event classes 
(INEL and NSD) introduced so far, results are also pre- 
sented for inelastic events with at least one charged parti- 
cle produced in the region \r)\ < 0.5, labeled as INEL>0. 
These values were obtained as the mean values of the cor- 
responding corrected multiplicity distributions for Neb. > 
(see Fig. H}. 

The consistency between data and model calculations 
varies with event class and the collision energy. PYTHIA 
tunes D6T and Perugia-0 significantly underestimate the 
charged-particle density in all event classes and at both 
energies. ATLAS-CSC tune, PHOJET, and QGSM are 
closer to the data and describe the average multiplicity 
reasonably well, at least for some of the classes and ener- 
gies listed in the Table [3] However, the relative increase in 
charged-particle density is underestimated by all models 
and tunes, most significantly for the event class with at 
least one charged particle in the central region (INEL>0). 
The increase predicted by all PYTHIA tunes is 16 % (14 % 
for PHOJET), whereas the observed increase is substan- 
tially larger (22.6 ± 0.7 ± 1.0) %. 

Figure [7] shows the centre-of-mass energy dependence 
of the pseudorapidity density in the central region. The 
data points are obtained in the \rj\ < 0.5 range from this 
experiment and from [40H48] . When necessary, corrections 
were applied for differences in pseudorapidity ranges, fit- 
ting the pseudorapidity distributions around rj = 0. 

Using parameterizations obtained by fitting a power- 
law dependence on the centre-of-mass energy, extrapola- 
tions to the centre-of-mass energy of yfs = 7 TeV give 
dN ch /dr] = 4.7 and dN ch /dr] = 5.4 for INEL and for NSD 
interactions, respectively. At the nominal LHC energy of 
y/s — 14 TeV, the same extrapolations yields dN^/dr/ = 
5.4 and dN ch /dr] = 6.2 for INEL and for NSD collisions, 
respectively. 



The multiplicity distributions of charged particles were 
measured in three pseudorapidity intervals at both ener- 
gies. These distributions, corrected as described above, are 
shown in Fig. [8] (left) and Fig. [8] (right) respectively, for 
y/s = 0.9 TeV and yfs = 2.36 TeV for NSD events. The 
difference between the multiplicity distributions for NSD 
and for INEL events only becomes significant at low mul- 
tiplicities (see Fig. IH]), as expected. 

In the two larger pseudorapidity intervals, small wavy 
fluctuations are seen at multiplicities above 25. While vi- 
sually they may appear to be significant, one should note 
that the errors in the deconvoluted distribution are corre- 
lated over a range comparable to the multiplicity resolu- 
tion (see Fig.H]). We studied the significance of these fluc- 
tuations assuming an exponential shape of the corrected 
distribution in the corresponding multiplicity range. Ap- 
plying the response matrix to this smooth distribution 
and comparing with the measured raw distribution, we 
find differences of up to two standard deviations in some 
of the corresponding raw data bins. Therefore, we con- 
clude that while the structures are related to fluctuations 
in the raw data, they are not significant, and that the un- 
certainty bands should be seen as one-standard-deviation 
envelopes of the deconvoluted distributions. Similar obser- 
vations for a different deconvolution method were made by 
UA5 in Q]5]. 

The multiplicity distributions were fitted with a Nega- 
tive-Binomial Distribution (NBD) and at both energies 
satisfactory descriptions were obtained, as shown in Fig. [5] 
Fitting the spectra with the sum of two NBDs, as sug- 
gested in , did not significantly improve the description 
of the data. 

A comparison of the data to the multiplicity distri- 
butions obtained with the event generators is shown in 
Fig.[inifor \rj\ < 1.0. At low multiplicities (< 20) discrep- 
ancies are observed at both energies and for all models. At 
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Fig. 9. Expanded views of the low-multiplicity region of corrected multiplicity distributions for INEL and NSD events, left for 
0.9 TeV and right for 2.36 TeV data. The gray bands indicate the systematic uncertainty. Distribution for NSD events are not 
normalized to unity but scaled down in such a way that the distributions for INEL and NSD events match at high multiplicities, 
which makes the difference at low multiplicity clearly visible. Left: data at y/s — 0.9 TeV. Right: data at y^s = 2.36 TeV. Note 
that for \r)\ < 1.0 and \rj\ < 1.3 the distributions have been scaled for clarity by the factor indicated. 
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Fig. 10. Comparison of measured multiplicity distributions for INEL events to models for the pseudorapidity range [77I < 1.0. 
Predictions are shown based on the PHOJET model [12] (solid line) and PYTHIA tunes: D6T [9] (dashed line), ATLAS-CSC [10] 
(dotted line), and Perugia-0 |11] (dash-dotted line). The error bars for data points represent statistical uncertainties, the shaded 
areas represent systematic uncertainties. Left: data at 0.9 TeV. Right: data at 2.36 TeV. For both cases the ratios between the 
measured values and model calculations are shown in the lower part with the same convention. The shaded areas represent the 
combined statistical and systematic uncertainties. 



high multiplicities and for the 0.9 TeV sample, the PHO- 
JET model agrees well with the data. The PYTHIA tunes 
D6T and Perugia-0 underestimate the data at high mul- 
tiplicities and the ATLAS-CSC tune is above the data in 
this region. At 2.36 TeV, ATLAS-CSC tune of PYTHIA 
and, to some extent, PHOJET are close to the data. The 
ratios of data over Monte Carlo calculations are very sim- 
ilar in all three pseudorapidity ranges and suggests that 



the stronger rise with energy seen in the charged-particle 
density is, at least partly, due to a larger fraction of high- 
multiplicity events. 

From these multiplicity distributions we have calcu- 
lated the mean multiplicity and first reduced moments 

C q = (N« h )/(N ch r, (5) 
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Table 4. Mean multiplicity and Cj-moments (J5]l of the multiplicity distributions measured by UA5 [19] in proton-antiproton 
collisions at y/s = 0.9 TeV, and by ALICE at yfs = 0.9 TeV and 2.36 TeV, for NSD events in three different pseudorapidity 
intervals. The first error is statistical and the second systematic. 





UA5 pp 


ALICE pp 




v^ = 


0.9 TeV 


y/s = 2.36 TeV 


\r)\ < 0.5 


(iV ch ) 


3.61 ±0.04 ±0.12 


3.60 ±0.02 ±0.11 


4.47 ±0.03 ±0.10 


c 2 


1.94 ±0.02 ±0.04 


1.96 ±0.01 ±0.06 


2.02 ± 0.01 ± 0.04 


c 3 


5.4 ±0.2 ±0.3 


5.35 ±0.06 ±0.31 


5.76 ± 0.09 ± 0.26 


Ct 


19 ± 1 ± 1 


18.3 ± 0.4 ± 1.6 


20.6 ±0.6 ±1.4 


H < i-o 


(iVch) 


7.38 ±0.08 ±0.27 


7.38 ±0.03 ±0.17 


9.08 ± 0.06 ± 0.29 


c 2 


1.75 ±0.02 ±0.04 


1.77 ±0.01 ±0.04 


1.84 ±0.01 ±0.06 


c 3 


4.4 ± 0.1 ± 0.1 


4.25 ±0.03 ±0.20 


4.65 ± 0.06 ± 0.30 


c 4 


14.1 ±0.9 ±1.2 


12.6 ±0.1 ±0.9 


14.3 ±0.3 ±1.4 


\V\ < 1-3 


(N ch ) 




9.73 ±0.12 ±0.19 


11.86 ±0.22 ±0.45 


c 2 




1.70 ±0.02 ±0.03 


1.79 ±0.03 ±0.07 


c 3 




3.91 ±0.10 ±0.15 


4.35 ±0.16 ±0.33 


c 4 




10.9 ±0.4 ±0.6 


12.8 ±0.7 ±1.5 



u 

£20 

i 15 

o 

% 5 



20 




15 
10 
5 




0.5 



to 



< 1.0 



UA5 



• Co 



» 



■■■□■— 

o : 



10 2 



10 J 



Energy s[s (GeV) 



Fig. 11. Energy dependence of the Cj-moments (J5J of the 
multiplicity distributions measured by UA5 QIJ] and ALICE at 
both energies for NSD events in two different pseudorapidity 
intervals. The error bars represent the combined statistical and 
systematic uncertainties. The data at 0.9 TeV are displaced 
horizontally for visibility. 



summarized in Table SJ For |^| < 0.5 and < 1.0 our 
results are compared to the UA5 measurement for pp colli- 
sions at y/s = 0.9 TeV [T!5]. Note that the mean multiplici- 
ties quoted in this table are those calculated from the mul- 
tiplicity distributions and are therefore slightly different 
from the values given in Table [3] The value of the pseudo- 
rapidity density obtained when averaging the multiplicity 
distribution for \rj\ < 0.5 is consistent with the value ob- 
tained in the pseudorapidity-density analysis. This is an 
important consistency check, since the correction methods 



□ UA5 0.2 TeV 

o UA5 0.9 TeV 

a ALICE 0.9 TeV 

o ALICE 2.36 TeV 




Fig. 12. Comparison of multiplicity distributions in KNO vari- 
ables measured by UA5 18,19] in proton-antiproton collisions 
at s/s = 0.2 TeV and 0.9 TeV, and by ALICE at s/s = 0.9 TeV 
and 2.36 TeV, for NSD events in \r)\ < 0.5. In the lower part the 
ratio between ALICE measurements at 0.9 TeV and 2.36 TeV 
is shown. The error bars represent the combined statistical and 
systematic uncertainties. 



in the pseudorapidity-density and multiplicity-distribution 
analyses are different. 

Our data are consistent with UA5 proton-antiproton 
measurements at 900 GeV (Fig. [5^ and Table E]). The en- 
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ergy dependence of the reduced moments C q , shown in 
Fig-El indicates a slight increase, which is not significant 
given the size of our systematic uncertainties. Systematic 
uncertainties are assumed to be uncorrelated between en- 
ergies. A similar conclusion about the shape evolution 
of multiplicity distributions can be drawn from Fig. 1121 
where we compare our measurements, plotted in terms 
of KNO variables, at the two energies and UA5 pp data 
at s/s = 0.2 and 0.9 TeV, for NSD collisions and pseu- 
dorapidity interval \r)\ < 0.5. While KNO scaling gives a 
reasonable description of the data from 0.2 to 2.36 TeV, 
the ratio between the 0.9 TeV and 2.36 TeV data shows a 
slight departure from unity above z — 4. 

7 Conclusion 

We report high-statistics measurements of the charged- 
primary particle pseudorapidity density and multiplicity 
distributions in proton-proton collisions at ccntre-of-mass 
energies of 0.9 TeV and 2.36 TeV with the ALICE de- 
tector. The results at 0.9 TeV are consistent with UA5 
pp measurements at the same energy. At both energies, 
our data are consistent with the CMS measurement, and 
compared to various models for which they provide further 
constraints. None of the investigated models and tunes de- 
scribes the average multiplicities and the multiplicity dis- 
tributions well. In particular, they underestimate the in- 
crease in the average multiplicity seen in the data between 
0.9 TeV and 2.36 TeV. At 0.9 TeV, the high-multiplicity 
tail of the distributions is best described by the PHOJET 
model, while at 2.36 TeV, PYTHIA tune ATLAS-CSC is 
closest to the data. 

The multiplicity distributions at both energies and in 
pseudorapidity ranges up to < 1.3 are described well 
with negative binomial distributions. The shape evolution 
of the multiplicity distributions with energy was studied 
in terms of KNO-scaling variables, and by extracting re- 
duced moments of the distributions. A slight, but only 
marginally significant evolution in the shape is visible in 
the data for z > 4, possibly indicating an increasing frac- 
tion of events with the highest multiplicity. This issue will 
be studied further using the data collected from forthcom- 
ing higher-energy runs at the LHC. 

Acknowledgements 

The ALICE collaboration would like to thank all its engineers 
and technicians for their invaluable contributions to the con- 
struction of the experiment and the CERN accelerator teams 
for the outstanding performance of the LHC complex. 

The ALICE collaboration acknowledges the following fund- 
ing agencies for their support in building and running the AL- 
ICE detector: 

— Calouste Gulbenkian Foundation from Lisbon and Swiss 
Fonds Kidagan, Armenia; 

— Conselho Nacional de Desenvolvimento Cientffico e Tec- 
nologico (CNPq), Financiadora de Estudos e Projetos (FI- 
NEP), Fundacao de Amparo a Pesquisa do Estado de Sao 
Paulo (FAPESP); 



— National Natural Science Foundation of China (NSFC), the 
Chinese Ministry of Education (CMOE) and the Ministry 
of Science and Technology of China (MSTC); 

— Ministry of Education and Youth of the Czech Republic; 

— Danish Natural Science Research Council, the Carlsberg 
Foundation and the Danish National Research Foundation; 

— The European Research Council under the European Com- 
munity's Seventh Framework Programme; 

— Helsinki Institute of Physics and the Academy of Finland; 

— French CNRS-IN2P3, the 'Region Pays de Loire', 'Region 
Alsace', 'Region Auvergne' and CEA, France; 

— German BMBF and the Helmholtz Association; 

— Hungarian OTKA and National Office for Research and 
Technology (NKTH); 

— Department of Atomic Energy and Department of Science 
and Technology of the Government of India; 

— Istituto Nazionale di Fisica Nucleare (INFN) of Italy; 

— MEXT Grant-in- Aid for Specially Promoted Research, Ja- 
pan; 

— Joint Institute for Nuclear Research, Dubna; 

— Korea Foundation for International Cooperation of Science 
and Technology (KICOS); 

— CONACYT, DGAPA, Mexico, ALFA-EC and the HELEN 
Program (High-Energy physics Latin- American-European 
Network) ; 

— Stichting voor Fundamenteel Onderzoek der Materie (FOM) 
and the Nederlandse Organisatie voor Wetenschappelijk 
Onderzoek (NWO), Netherlands; 

— Research Council of Norway (NFR); 

— Polish Ministry of Science and Higher Education; 

— National Authority for Scientific Research - NASR (Auton- 
tatea Nationala pentru Cercetare Stiintifica - ANCS); 

— Federal Agency of Science of the Ministry of Education and 
Science of Russian Federation, International Science and 
Technology Center, Russian Academy of Sciences, Russian 
Federal Agency of Atomic Energy, Russian Federal Agency 
for Science and Innovations and CERN-INTAS; 

— Ministry of Education of Slovakia; 

— CIEMAT, EELA, Ministerio de Education y Ciencia of 
Spain, Xunta de Galicia (Conselleria de Education), CEA- 
DEN, Cubaenergia, Cuba, and IAEA (International Atomic 
Energy Agency); 

— Swedish Reseach Council (VR) and Knut & Alice Wallen- 
berg Foundation (K AW) ; 

— Ukraine Ministry of Education and Science; 

— United Kingdom Science and Technology Facilities Council 
(STFC); 

— The United States Department of Energy, the United States 
National Science Foundation, the State of Texas, and the 
State of Ohio. 



References 

1. ALICE Collaboration, K. Aamodt et al., J. Instrum 3, 
S08002 (2008) 

2. ALICE Collaboration, K. Aamodt et al., Eur. Phys. J. C 
65, 111 (2010) 

3. L. Evans, P. Bryant (editors), J. Instrum. 3, S08001 (2008) 

4. CMS Collaboration, V. Khachatryan et al., JHEP 2010 
02 041 (2010) 



19 



5. ATLAS Collaboration, G. Aad et al, submitted to Phys. 
Lett. B, CERN preprint CERN-PH-EP/2010-004 (2010), 
arXiv:1003.3124 [hep-ex] (2010) 

6. A.B. Kaidalov and M.G. Poghosyan, to be published in 
Eur. Phys. J. C, arXiv:0910.2050[hep-ph] (2009) 

7. T. Sjostrand, Comput. Phys. Commun. 82, 74 (1994) 

8. T. Sjostrand, S. Mrenna, P. Skands, J. High Energy Phys. 
2006, 05 026 (2006) 

9. M.G. Albrow et al. (Tev4LHC QCD Working Group), 
arXiv:hep-ph/0610012 (2006), D6T (109) tune 

10. A. Moraes (ATLAS Collaboration), ATLAS Note ATL- 
COM-PHYS-2009-119 (2009), ATLAS CSC (306) tune 

11. P.Z. Skands, Multi-Parton Interaction Workshop, Perugia, 
Italy, 28-31 Oct. 2008, arXiv:0905.3418[hep-ph] (2009), 
Perugia-0 (320) tune 

12. R. Engel, J. Ranft, S. Roesler, Phys. Rev. D 52, 1459 
(1995) 

13. Z. Koba, H.B. Nielsen and P. Olesen, Nucl. Phys. B 40, 
317 (1972) 

14. A. Breakstone et al., Phys. Rev. D 30, 528 (1984) 

15. UA5 Collaboration, G.J. Alner et al., Phys. Lett. B 167, 
476 (1986) 

16. UA5 Collaboration, G.J. Alner et al., Phys. Lett. B 138, 
304 (1984) 

17. UA1 Collaboration, C. Albajar et al., Nucl. Phys. B 335, 
261 (1990) 

18. UA5 Collaboration, G.J. Alner et al., Phys. Lett. B 160, 
193 (1985) 

19. UA5 Collaboration, R.E. Ansorge et al., Z. Phys. C 43, 
357 (1989) 

20. CDF Collaboration, F. Rimondi et al., in 23rd Int. Symp. 
on Multiparticle Dynamics (1993), Aspen, Colorado, 12-17 
Sep. 1993, Word Scientific, 400, (1994) 

21. A. Giovannini and R. Ugoccioni, Phys. Rev. D 59 094020 
(1999), erratum Phys. Rev. D 69, 059903 (2004) 

22. A. Giovannini and R. Ugoccioni, Phys. Rev. D 60, 074027 
(1999) 

23. J.F. Grosse-Oetringhaus and K. Reygers, to be published 
in J. Phys. G, arXiv:0912.0023[hep-ex] (2009) 

24. ALICE Collaboration, K. Aamodt et al., JINST 5, (2010) 
P03003 

25. ALICE Collaboration, R. Brun et al., Nucl. Instrum. Meth. 
A 502, 339 (2003) 

26. R. Brun et al., 1985 GEANT3 User Guide, CERN Data 
Handling Division DD/EE/841 and 1994 CERN Program 
Library Long Write-up, W5013, GEANT Detector De- 
scription and Simulation Tool 

27. UA5 Collaboration, R.E. Ansorge et al., Z. Phys. C 33, 
175 (1986) 

28. UA4 Collaboration, D. Bernard et al., Phys. Lett. B 186 
(1987) 227 

29. A.B. Kaidalov and M.G. Poghosyan, arXiv:0909.5156 [hep- 
ph] (2009) 

30. N.A. Amos et al, Phys. Lett. B 301, 313 (1993) 

31. F. Abe et al., Phys. Rev. D 50, 5535 (1994) 

32. S. Klimenko, J. Konigsberg and T.M. Liss, Fermilab 
preprint, Fermilab-FN-0741 (2003) 

33. CDF Collaboration, T. Affolder et al., Phys. Rev. Lett. 
87, 141802 (2001) 

34. PHOBOS Collaboration, B.B. Back et al., Phys. Rev. Lett. 
85, 3100 (2000) 

35. V. Blobel, in 8th CERN School of Comp. - CSC'84, 
Aiguablava, Spain, 9-22 Sep. 1984, CERN-85-09, 88, 
(1985) 



36. G. D'Agostini, Nucl. Instrum. Meth. A 362, 487 (1995) 

37. G. D'Agostini, CERN Report CERN-99-03 (1999) 

38. J.F. Grosse-Oetringhaus, PhD thesis, University of Minis- 
ter, Germany, CERN-THESIS-2009-033 (2009) 

39. G. Cowan, in Advanced Statistical Techniques in Parti- 
cle Physics, Durham, England, 18-22 Mar 2002, Durham 
Univ., 248 (2002) 

40. UA5 Collaboration, G.J. Alner et al., Z. Phys. C 33, (1986) 
1 

41. W. Thome et al., Nucl. Phys. B 129, 365 (1977) 

42. UA5 Collaboration, K. Alpgard et al., Phys. Lett. B 112, 
183 (1982) 

43. M. Ambrosio et al., AIP Conference Proceedings 85, 602 
(1982) 

44. PHOBOS Collaboration, R. Noucier et al., J. Phys. G 30, 
S1133 (2004) 

45. STAR Collaboration, B.I. Abelev et al., Phys. Rev. C 79, 
034909 (2009) 

46. UA5 Collaboration, G.J. Alner et al., Phys. Rep. 154, 247 
(1987) 

47. UA1 Collaboration, C. Albajar et al., Nucl. Phys. B 335, 
261 (1990) 

48. CDF Collaboration, F. Abe et al., Phys. Rev. D 41, 2330, 
(1990) 

49. Ch. Fuglesang, in Multiparticle Dynamics - Festschrift for 
Leon Van Hove, La Thuile, Italy, 10-22 Mar. 1989, World 
Scientific, 193 (1989) 



